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Abstract: Following Green’s function technique and equation of motion 
method, the coexistence of superconductivity (SC) and itinerant ferromagnetism 
(FM) is investigated in a single band homogenous system. Self consistent 
equations for SC and FM order parameters, A and m or I respectively are 
derived. It is shown that there generally exists a coexistent (A ^ 0, and m 
or I ^ 0) solutions to the coupled equations of the order parameter in the, 
temperature range 0 < T < min (T c T FM ), where T and T rM are respectively the 
superconducting and ferromagnetic transition temperatures. Expressions for 
specific heat, density of states, free energy and critical field are derived. The 
specific heat has linear temperature dependence as opposed to the exponential 
decrease in the BCS theory. The density of states for a finite m increases as 
opposed to that of a ferromagnetic metal. Free energy study reveals that FM- 
SC state has lowest energy than the normal FM state and therefore realized 
at low enough temperature .Effect of small external field is also studied. The 
theory is applied to explain the observations in uranium based intermetallics 
systems UCoGe and Ulr. The agreement between theory and experiments is 
quite encouraging. 

Keywords: Itinerant ferromagnetism, Green’s function, Superconductivity, 
Energy spectra and density of states, BCS Hamiltonian, 

Hubbard Hamiltonian. 


Journal of Nuclear 
Physics, Material 
Sciences, Radiation 
and Applications 
Vol. 2, No. 1 
August 2014 
pp. 33-71 


33 


CHITKARA 

UNIVERSITY 





1. INTRODUCTION 


Kakani, S. 
Nuwal, A. 
Kakani, S.L. 


The recent discovery the coexistence of superconductivity (SC) and itinerant 
ferromagnetism (FM) in Uranium based intermetallics compounds: UGe 2 
[1], URhGe [2], Ulr [3] and UCoGe [4] have attracted a great interest and 
revived the interest in the old problem of coexistence of superconductivity and 
ferromagnetism [5-9]. In these Uranium intermetallics magnetism has a strong 
itinerant character and both ordering phenomena are carried by the same 5f 
electrons. The experiments indicate that in these materials: 

(i) The superconducting phase is completely covered within the ferromagnetic 
phase and disappears in the paramagnetic phase. 

(ii) The ferromagnetic order is stable within the superconducting phase. 

(iii) The specific heat anomaly associated with the superconductivity in these 
materials appears to be absent. The specific heat depends on the temperature 
linearly at low temperature [10]. 

Numerous theoretical studies have tried to answer the question what the 
attractive forces are between the electrons leading to the formation of Cooper 
pairs [11-32], The general features of the proposed models are: 

(i) The superconducting pairing of the conduction electrons is mediated by 
spin fluctuations rather than by phonons, as is the case in conventional 
superconductors. 

(ii) In the superconducting state quasiparticles form Cooper pairs in which the 
spins are parallel (S = 1) in contrast to conventional superconductors with 
opposite spin (S = 0). 

(iii) The ferromagnetism is itinerant and therefore carried by the conduction 
electrons. This arises from a splitting of the spin-up and spin-down band. 
A consequence is that the ferromagnetism and superconductivity is carried 
by same electrons. 

Early theories proposed magnetically mediated spin - triplet 
superconductivity [27]. Suhl [19], Zhou and Gong [24 ], Mineev and 
Champel [17] have also favoured magnetically mediated spin - triplet 
superconductivity in these systems. Following Suhl’s work [19], Abrikosov 
[20] have shown that this mechanism can lead only to an s - wave order 
parameter. Superconductivity appears together with ferromagnetism but 
persists only until the ferromagnetism is weak. Watanable and Miyake [13] 
are of the opinion that the superconductivity in these systems is mediated 
by the coupled charge density wave (CDW) and spin density wave (SDW) 
fluctuations originating from the CDW ordering of the majority - spin band. 


34 





Powell et al. [25] have considered spin singlet and spin triplet states with 
either equal spin pairing (ESP) or opposite spin pairing (OSP) states. They 
found that gap equations for the singlet state reproduced the Clogston - 
Chandrasekhar limiting behaviour [33,34] and the phase diagram of the 
Baltensperger - Sarma equation [35,36]. They also showed that the singlet 
gap equation leads to the result that the superconducting order parameter is 
independent of exchange splitting at zero temperature. They further found 
that OSP triplet states showed a very similar behaviour to the singlet state 
in the presence of exchange splitting. Blagoev et al. [26] have studied the 
problem in terms of diagrammatic many - body theory, on the basis of the 
itinerant electron model. Machida and Ohmi [11] have given arguments for 
the pairing to be of p - wave, normally spin - triplet nature. It is not obvious 
from these considerations why the superconductivity is observed only in 
the ferromagnetic phase. To explain this Sandeman et al. [14 ] proposed 
a density of states effect that exists only in the ferromagnetic phase, as 
the source of the superconductivity in UGe v Kirkpatric et al. [15,16] 
proposed an explanation for the observed phase diagram that is based on an 
enhancement of the longitudional spin susceptibility in the ferromagnetic 
phase by magnons or magnetic Goldstone modes. Most theories which 
describe ferromagnet close to a quantum phase transition have predicated 
that the superconducting transition temperature, T c , should be at least as 
high in the paramagnetic state as it is in the ferromagnetic state. These 
theories have considered an electronically three - dimensional ferromagnet, 
either magnetically isotropic [21] or uniaxial [23]. 

Kirkpatrick et al [15,16] have predicated an enhancement of the 
superconducting transition temperature T c in the ferromagnetic regime 
from the coupling of magnons to the longitudinal magnetic susceptibility. 
However, the ferromagnetic state of these systems is highly anisotropic - at 
4.2 K and an external magnetic field of 4T, the easy - axis magnetization 
is 20 to 30 times that along either of the other crystallographic axes [37] - 
so transverse modes seem unlikely to explain the exclusively ferromagnetic 
superconductivity in these systems. 

Gorski et al. [3 8] have studied the possibility of simultaneous coexistence of 
superconductivity and ferromagnetism within the framework of the extended 
Hubbard model. The main driving forces for these phenomena considered 
by these authors are kinetic interactions. They reported that even weak 
feiTomagnetism, if generated by the band shift, destroys the superconductivity. 
The ferromagnetism created by a change of bandwidth can coexist with the 
singlet superconductivity. However, it seems difficult to understand the 
simultaneous coexistence of superconductivity and ferromagnetism in these 
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systems in terms of these theories and require a novel concept to understand : 
(i) how the same band electrons are responsible for both the superconductivity 
and ferromagnetism, and (ii) why the superconductivity occurs only in the 
ferromagnetic phase, i.e. what is the microscopic relationship between these 
two antagonistic phenomena? 

The theories of itinerant ferromagnetic superconductors stated above 
provide a basic understanding of the new phenomenon in these systems, but 
a microscopic theory and a satisfactory mechanism seems still lacking. 

In the present paper, we propose a simple microscopic theory of 
ferromagnetic superconductors which explain satisfactorily the observed 
features of these systems and also predict the qualitative behaviour of energy 
spectra and density of states function N(u), superconducting order parameter 
(A), magnetic order parameter (m), electronic specific heat, free energy , 
critical field and effect of small external field. In this class of materials, Ulr 
holds a special place because it is a system in which the crystal structure 
lacks inversion symmetry. 


2. THE MODEL HAMILTONIAN 

We propose the following model Hamiltonian: 


with 


where 


H — H 0 

Hbcs + H h 

(1) 


(K)-»} C + Ka C Ka 

(2) 

K,a 




{<e(k)-h b 1vj} 

(3) 


1 for the spin up 

a = \ 

[-1 for the spin down 

h (= // ( . H) is small external field in tesla. C + Kn (C Ka ) are the fermion creation 
(annihilation) operators for the Block states, | K, a >, where K is the conduction 
electron wave vector and a is its spin. G(K) is the single electron energy in 
the f-band and measured from the bottom of the band and // is the chemical 
potential. The operators C^,C Ki7 satisfy the usual fermion anti - commutation 
rules. 


Hbcs — ^^(C_kjC kt +C kt C_ K | ) (4) 

K 
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is the BCS Hamiltonian [39]. A is superconducting order parameter and 
is chosen real and positive. A is only a weak function of K and therefore we 
have neglected its K dependence. A is determined self - consistently from 
the gap equation 

A “ CrtQkj. > (5) 

-O K 

|g| (|g| > 0) is the phonon - mediated electron - electron coupling constant 
having the dimensions of energy. The summation over K is limited by the 

wave vector corresponding to Debye’s energy hoj D at the Fermi surface. <- 

> indicates the thermodynamic average. 
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H h = 


u 

N 


E r+ r+r r 

'-'ATT Kl Kj'-KT 

K, K’ 


( 6 ) 


is the Hubbard Hamiltonian [40]. U is the repulsive Coulomb interaction, 
which is responsible for the ferromagnetic transition. 

Combining equations (2), (4) and (6), the model Hamiltonian (1) can be 
expressed in the following mathematical form. 


H = -M-AV^}CL C K.- A E]( C -n C n + C+ k£ + -ki) 


U 


+ T/ : K|C K -|C K X 


We define 


K| 


11 a ~ ivr X/ < > 

iN K 


and use the following notations. 

n = n | +n j= (C|C T +C|C 1 


(7) 

( 8 ) 

(9) 


m = nt-n|=(c;C T -C;C i ) 


( 10 ) 


n mcr _ , , ,.,. 

n .=- + — (^ = ±D (ID 

Where n is the total number of electrons per unit cell, /i fi is Bohr magneton 
and p is magnetic permeability of free space. It is proposed to solve the 
Hamiltonian (7) with the help of Green’s function technique using the equation 
of motion method. 
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2.1 GREEN’S FUNCTIONS 

In order to study the physical properties of itinerant ferromagnetic 
superconductors, we define electron Green’s functions for the conduction band 
electrons which are subjected to the ferromagnetic and superconducting instability. 
The electron Green’s functions are defined in the usual way as follows: 

G n (K”, K’t) = «C K ,,;C+ t » u 

= -ifl(t) < [C K „ t (t),C+ T » (12) 

G u (K”,K’,t) = «C K „ r ;C+ i » u 

= -iO(t) < [C K , u (t),C^(0)]> (13) 


F+(K",K’,t) = <<C ;C+ T » W 


- K"J 


= 10 (t) < C + (t);C+ T (0) 


-K"t 


> 


(14) 


F-n (K",K’,t) = «C ;C K1 » W 


- K"t 

i 0 ft) < 


C K .. T (t); C,,(0) 


> (15) 


Following equation of motion method, we obtain following expressions for 
Green’s functions, 


«c Kt ,c; T »= 


5 kx.[ w + ( G k -ti-M B ha)-UnT] 


[{w ^ ( e K - / i B hcr ) - Un 1} {uj + (€ k -/x - n B ha) - Un |} - A 2 ] 

(16) 


«c k .,c k ..», 


5 KK.[ w + ( e K -F - M B h <0 + Un |] 


[{a) + (e K -fi - Mb ho-) - Un t}{w + (e K -A( - M B ho-) + Un f}-A ; ] 


(17) 


+ + -AlS 

«c ,C >> =-—- 

[{<^ — (e K l}{w + (e K — /qha) — Un T} A” j 


(18) 
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2.2 CORRELATION FUNCTIONS 

Using the relation [5,7,41] 


oo «A(t);B(t’) -«A(t);B(t» 

<B(t’)A(t)>=lim— f - 

W W 2tt e BtJ -l (19) 


"OO 


x[exp[iw(t — t’)}]dw] 


and employing the identity, 

lim 1 

>0 


u) -\- i G —E„ u} — i G —E„ 


= —27r i 8 (cu-E K ) 


( 20 ) 


one obtains the correlation function for the Green’s functions given by 
equations (16) to (18) as 


<x 


< C+ C >= f(a ) + —L {f(a ) _ f(a )} 
iv | iv | 2 a, — au 1 2 


1 2 


+ 


( e K + 


a l a 2 


{f( ttl )-f(a 2 ) 


or 


<C+ .C^ . > = f(aj- 


«1 + —// - /ighcr j + Ura l 


'K| K| 


a r a 2 


( 21 ) 


x{f( ai )-f(a 2 )} 


Where 


And 


Um 


Un 


^ a r^ 2 +1 i e K ~ + — 


Um 


a, = 


~ + j( e /r ~ 


( 22 ) 


and /(aj) and /(a,) are Fermi functions. Similarly, one obtains the other 
correlation functions as 


< C + C > = 


a, •+ a, ’+ (g K - ,.„H + Un l . (<> . } _ f} j (2J) 

rv ’— rv ’ t J 
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where a { = ar >d a 2 = — a, and f( a { ) and f (a 2 ) are Fermi functions. 

< C;' K , T C^ >= - —[{/(a,) - f(at 2 )}] (24) 

Using these Green’s functions and correlation functions the expressions for 
various physical properties of itinerant ferromagnetic superconductors can be 
derived. 


3. PHYSICAL PROPERTIES 


3.1 SUPERCONDUCTING ORDER PARAMETER (A) 

A is determined self - consistently from the gap equation 

^ = "TrXY C K jC-K| > (25) 

™ K 

|g| (|g| > 0) is the phonon - mediated electron - electron coupling constant 
having the dimensions of energy. The summation over K is limited by the 
wave vector corresponding to Debye’s energy Huj d at the Fermi surface. < 

-> indicates the thermodynamic average. By substituting the correlation 

function given by equation (24). 


or 


— = I 

N(o) J i 

—=r 

NY rO J M 


d€. 


1 

1 

(a, — a 2 ) 

+1 


1 


• /i+/ku D 

. de K 

|g|N(o) J '‘-'"d 


1 


4 '\| + _/z " Mb/kj ) + y} 


Ba, 8a, 

tanh -—- — tanh -—- 



(26) 


Where a and a, are given by equation (22). 

Putting m = 0 and U = 0 in equation (26), we observe that it reduces to the 
standard BCS equation [39], 
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1 r d dx , x 

— -= I —tanh- (27) 

|g|N(o) J o x 2 k B T c 

where k B is Boltzmann constant and T c is superconducting critical transition 
temperature. 

Analyzing the expression (26) , we see that the presence of magnetization 
reduces the volume of the phase space that is available for the Cooper pair. 
This leads to the decrease in gap at T = 0 and hence T. will also be reduced 
compared to the BCS value for given values of other relevant parameters. 
Obviously, the increase in magnetization pushes the gap away from the 
Fermi surface and reduces its magnitude also [42] 

Equation (26) is the general expression for the superconducting order 
parameter A for itinerant ferromagnetic superconductors. 
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3.2 MAGNETIZATION (m) 

From equation (10), we have - 

m=n|-n j 


Obviously, the fermions which form Cooper pairs are the same as those 
responsible for spontaneous magnetization. 

Using equation (8), the above relation takes the form. 

"■=^E<Q,c K1 >-4E<ct,CK,> (2» 

Using equations (21 and (23, the equation (28 takes the form. 


where 


m =-£ 

NY 


1 

1 

[e^ E -+l 

e^+ll 


E = E„ 


E+ = E k+ I 

are the quasipartile energy relations. Here, 

Urn 


Isl¬ 


and 



Un 

| A + 

(£ K -n - n B ha) + — 


(29) 

(30) 

(31) 

(32) 

(33) 
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Now, one can write the magnetic energy splitting I between the spin - up 
and spin - down band as 


2 2N k 


(34) 


or 


1= “ ^[f(E-)-f(E + )] 

, 4- K 


K 


(35) 


where u = and f (E) is Fermi function. 

Equation (26) is the BCS gap equation in the presence of a finite magnetization 
m. Equation (35) is the equation which determines the magnetization when 
the superconducting state may be present. Here the superconducting gap 
A and the magnetic energy splitting I between the spin-up and spin-down 
bands are independent order parameters. The system of equations (26) and 
(34) determines the phase where the superconductivity and ferromagnetism 
coexist. 

The equations (26) and (34) have the following solutions : 

(i) normal or paramagnetic state A = 0,1 = 0; 

(ii) ferromagnetic state A = 0,1 ^ 0; 

(iii) superconducting state A ^ 0,1 = 0; and 

(iv) coexistent state A ^ 0,1 ^ 0. 

From equations (26) and (34), we find that magnetization m and superconducting 
order parameter A have dependencies such as m =m (g, U, T) and A = A (g, 
U, T). Moreover, they are coupled with each other through the distribution 
functions, which are also functions of m and A. From the numerical solutions 
of these two coupled equations (26) and (34), we will get the values of m and 
A to be used for the calculation of the physical properties. 

The ferromagnetic transition temperature T FM can be obtained from 
Equation (29). Equation (29) cannot be solved exactly and therefore we 
evaluate it numerically. 

3.3. SPECIFIC HEAT (C s ) 

The electronic specific heat per atom of a superconductor is determined from 
the following relation [6,8,43]. 


d 1 

— At ~ Ms^ 17 ) < C K jC K | > 


(36) 
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Changing the summation over K into integration and substituting the correlation 
function < C^ T C K1 > , one obtains, 
or 


C! 

T 


2N(o) 1 /3a 2 exp(/3a 2 ) 

N T 2 Jo k [exp(/?a 2 ) + l] 2 


+ 


/?(«!+ E K —/r B hcr + Un |)(e k ~^crh) 


2 + 1 ( e * + 


Un 


x 


cq exp(/3aj) 
[expC/?^) +1] 2 


a 2 exp(/5a 2 ) 
[exp(/3a 2 ) + l] 2 


(37) 
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One can solve equation (37) numerically for electronic specific heat. 


3.4 ENERGY SPECTRA AND DENSITY OF STATES 

If we compare the quasiparticle energy spectra given by equations (30) and (31) 
with the BCS energy spectrum E = +A 2 , we note that the modification 

is due to the presence of the ferromagnetic energy. We can easily verify that 
the gap appears not at the Fermi level, rather it is pushed up. Now, the energy 
spectrum is symmetrical around E = Um, contrary to the BCS case where it 
was symmetric around E = 0. 

Now, we derive the expressions for the density of states. Density of states 
is one of the most important function to be determined and the one which is 
most susceptible to experimental verification as a function of the excitation 
energy E K [6-8, 39]. 

The density of states per atom per spin, N(u) is given by [44] 


N(lu) = lim ^ 1 y^G n [(/r,cu + / E) — G u (K,uj — i E)] 


(38) 


Here G (K, uj) is the one particle Green’s function given by equation 
(17). From equations (38) and (17), using following relation: 


lim 

e-*o 


1 

i 

cv + / G — E k 

co — i E — E k 


= 2n iS ( oj — E k 


(39) 


we obtain 
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N(u>) = N( 0) 


(u> + n B hc r) + (Um/2) 
((cj + ^ fi /7cr) + (?7m/2)j —A 2 


(40) 


N(uo) _ (w + ^ B ha) + I 


(41) 


((u; + + /) — A 2 ^ 


In the limit I —> 0, equation (41) reduces to 




forte > 0 


(42) 


N(0) 


0 otherwise 


which is same as described by the BCS theory [39]. 

The important effect of the magnetization is that it modifies the density 
of states, which affect the superconducting gap. We evaluate equation (42) 
numerically for obtaining density of states. 

3.5 FREE ENERGY 

Alongwith the nonzero solution, describing the coexistence of ferromagnetism 
and superconductivity, there is a solution of equations (26) and (35) with a 
vanishing gap describing normal ferromagnetic state. For the transition from 
the normal ferromagnetic state to the superconducting - ferromagnetic state to 
take place the energy of the former state must be lower than the energy of the 
latter state. One can calculate the difference between two energies using the 
relation [44-53], 



(43) 


lis form is particularly convenient because equation (43) expresses 
as a function of A (T), and direct substitution leads to 



V 




(45) 
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Putting U = 0 and h = 0 in equation (45), it reduces to standard BCS 
expression [11], 


Fsf-Fnf 

V 


n n A 

= N (°)/„ d ^/„ dA(A) ; 


d 

tanh 


dA 

E 


(46) 


Where E = ^A 2 +Q 
Writing equation (45) as 


F -F 

SF_NF 

V 


= N 


<») />./>*■> 


d 

dA 


2JA- + 


e K ~V B vh + - 


Un 

2 


(47) 


tanh ^ 


or 

p SF - f nf n (°) r ^ 


v 


0 d€ J o dA(A) 


dA 


-tanh 


/3a, 




Un 

2 


iV(0) 

r* hu D pA d 

(\ (^ \ HA(A ) 


1 

tanh 

2 “ 

0 j 0 dA 

L 

Un) 

2 2 



If + 

e K -P B h<J + — 

v. 2 





A 2 



1 

Um I , 

Un 

2 

tanh — /3 

+ Ja + 

G k + - 


2 

2 V 

2 
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1 

tanh — /3 

2 





( 48 ) 



A 



where the second and fourth line is obtained with partial integration. The first 
and third terms together is just similar to the right side of the BCS gap equation 
[ 11 ], 


1 _ f^ D dG * tnnh ^^ 2 ) 

g\N(0) Jo (4+A 2 ) 1/2 2 


(49) 


The second and fourth terms can be simplified by changing variables from 


A to t =, A 2 


i Un 

G K +fJ- B hcr + — 


. In this way, one obtains 


F sf- F np N(0)A 2 


V 


/ , ru: 

0 


de. 


, tanh 


A 2 


e K 


Un 


1 -I3a, 


P 


) n 1 „ 

f o de * £n 


cosh 


P 


Um 

2 


A 2 


G k 


Un 


cosh 


P 


Um 

2 


e K ' 


Un 
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N(0)A 2 r^ B 


/ , fu 
0 


d& 


A 2 + G K —/^ B h(j+ 


UnV 


- tanh 


I/Ja, 


4N(0) 


(3 


/ «“ D 
0 


d<5. 




cosh — 


■T-f’ 


, Un 
e K _ M B (Th + “T" 


cosh 


/? 


Um 

, Un 'll 

2 

£ k -p B erh + — 
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I r f » 
J o 


de. 


In 


1 + e 


Um 

2 


A2+ ( e K +M B ffh + I f) 


1 

Um I 

Un 

[Um Uni 

+ ~ P 

-— + Ja + 

£ k +- 

- f h e K —P B h(7 4-[ 

2 

2 \ 

l 2 J 

l 2 2 JJJ 


4N(0) pH 

+- de K 111 

2 J 0 K 


l + exp-j-^(--^21+e K -p B ah + -^- 


4N(0) rto. 


(3 


-f 

J o 


d In 


1 + exp 


~f3 




i' 


Um | , 

Win) 

2 

Um Un 


- —' ~J A + 

<E k -n B ha + — 

— ‘ 

—— e * ~^B ha + — 


2 \ 

2 


[2 2 J 



l p _ 

- I d€ In 

/3 J° 


-/? 


1 + e 


Um 


i Un 

g k^b CTh + T 


Since fiuj D >> k H T , the second and fourth integrals on the right may be 
extended to infinity. An easy calculation then shows that it equals the first 
temperature - dependent correction to the thermodynamic potential in the 
normal ferromagnetic state [11]. We have 
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4 N{o)V r°° 

Jo 


In 


-0 


1 + e 


Um 


G 


K 


i Un 

—fi ah H- 

B 2 


= ^N(0)VTT 2 (k B Tf 

«[Fnp(T)-F nf (o)] 

In addition, it is readily verified that [11] 


( 51 ) 


r 

Um 1 2 

Un ] 

2 


n d£ K 

+Ja + 

£ k + — 


- (e K —M B ch) 

JO 

2 V 

2 



‘ -N( 0) 


1 ., ., 2hcj n 

-A 2 +A 2 In- 2 

2 A 


(52) 


= - - A(o)A 2 - — - N( 0)A 2 In ^ 
2 g A 


After simplification, one obtains the expression for free energy difference between 
superconducting - ferromagnetic state and normal ferromagnetic state as 


-—= — — A(0)A 2 - 
V 2 

- 4N(0)k B T d € k 

+ ~TT 1 N(0)(k T) 2 
3 


riV(0)A 2 ln A ° 


A 


In 


1 + e 


Um 7 ( UnY 

+^A +^ K -fi B ha +—j 


(53) 


Where A 0 is the zero temperature gap, whose value can be determined from 

/3A 0 = ire~' = 1.7i (54) 

which is a universal constant independent of the particular material. 

On further simplification, equation (53) finally takes the form 
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+ -7r 2 A(0)(£ B T) 2 
3 


This shows that the superconducting ferromagnetic state has lower energy 
than the normal ferromagnetic state and therefore will be realized at low 
enough temperature [39,54 ]. , we evaluate numerically equation (55) for 
free energy. 


3.6 LOW TEMPERATURE CRITICAL FIELD (H ) 

The critical field H c (T) called the thermodynamic field for the transition 
between the normal ferromagnetic to superconducting ferromagnetic state is 
given by the difference of free energies [46]. 


He 

Sir 


= ~ F sf - F nf 


A combination with equation (55), equation (56) yield 


(56) 


H = (8 7r V) A 


- —iV(0)A" -tV(0)A' In 
2 


(4,1 

4iV(0) 7 

. A 

P 1 


d P In 


(57) 


1 + exp 


-P 


Um 


- + JA 2 + 


Un 

-vM + — 
2 


1 2 A(0) 


3 /f 


In the limit U —> 0 and h —> 0, equation (57) gives the low temperature critical 
field as 

2 ^ \ 2 ' 

H c (T) — [4 7r A (0) Ap ] " 


e 21 

l — 

3 


c 


(58) 


; H C (0) 


1-1.06 
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T, 
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T —> 0 
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H c (0) = [47riV(0)A^] 1/2 (59) 

is the critical field at T = 0. One can solve numerically equation (57) to get the 
behaviour of critical filed. 

4. NUMERICAL CALCULATIONS 

4.1 Superconducting and Magnetic order parameters 

Expression (26) is the expression from which we can study the variation 
of superconducting order parameter with temperature for ferromagnetic 
superconductors. Using the values of parameters given in tables 1 and 2 for 
UCoGe and Ulr and solving numerically, we obtain values of A at various 


Table 1: Values of Various Parameters for the System UCoGe 


Parameter 



Value 

References 

* Crystal Structure 



Orthorhombic TiNiSi structure 

68,69 




a = 6.845 A 
b = 4.206 A, 
c = 7.222 A 





Space group P 


* Curie temperature (T ) 



3 K 

4,59,70-71 

* Superconducting transition 


0.8 K 

4,59,70-71 

temperature (T) 





* Electronic Specific heat coefficient (T) 

0.057 J K- 2 mok 1 

72 

* Phonon energy ( hu> D ) 



1.7 x 10" 21 J 

9,73-75 

* Repulsive Coulomb energy U 


0.169 x 10- 19 J 

9,73-75 

* Number of electrons per unit cell (n) 


4 

68 

* Magnetization per atom 



= 0.7 x 10" 2 

74 

* Magnetic energy splitting 

' Urn 

2 


S 1.92 x 10 -21 J/atom 

74 

between spin up and spin down bands 

* BCS attractive interaction strength |g| 

3.8 x 10" 19 J/atom 

73-75 

* Density of States at the Fermi Surface, 

= 0.24 x 10 19 

73-75 

N(0) 



per Joule per atom 



Fermi Wave Vector, K 1.62 A -1 73-75 
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Table 2: Values of Various Parameters for the System Ulr 


Parameter 


Value 

References 

• Crystal structure 



Monoclinic PbBi type 
structure 

(space group P2 X ) 
a = 5.62 A, 
b= 10.59 A, 
c = 5.60 A, 

3 

• Superconducting Transition Temperature 
(T). 

0.14 K 

3 

• Curie Temperature (T^) (with low- 
temperature ordered moment) 


46 K at ambient pressure 

75,76 

• Electronic specific heat coefficient (T) 

0.049 J K- 2 mol" 1 

3 

• Phonon energy ( hlO D ) 



1.68 x 10- 21 J 

9,73-75 

• Repulsive Coulomb energy U 


0.158 x 10" 19 J 

9,73-75 

• Density of States at the Fermi Surface (N(0)) 

0.22 X 10 19 
per Joule per atom 

9,73-75 

• BCS attractive interaction strength |g| 

3.9 X 10" 19 J/atom 

9,73-75 

• Number of electrons per unit cell (n) 


8 

9,73-75 

• Magnetization per atom (m) 


= 0.71 x 10" 2 

74 

• Fermi wave vector, K 

• Magnetic energy splitting 

between spin up and spin c 

' Um 


1.62 A" 1 

1.68 x 10" 21 J/ atom 

9,73-75 

4,75,76 

2 

own banc 

s 
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temperatures below T c as shown in Table 3 and 4 and variation is shown in figs. 

1 and 2 respectively. From these graphs, we note that 

(i) Superconducting transition temperature for UCoGe and Ulr are 0.8 K and 
0.14 K respectively. 

(ii) A versus T variation is just similar to BCS. 

(iii) In comparison to BCS value of T c (T c BCS ), T c for ferromagnetic 
superconductors is slightly lower, i.e. magnetization suppresses 
superconductivity. This result is in agreement with Karchev [10,55], 
Karchev et al. [56] and Linder and Sudbo [54]. 
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Table 3: Superconducting order parameter (A) for UCoGe 


Temperature (T) 

K 

Superconducting order parameter 

A = x x 10" 21 J 

Theory 

h = 0.4 T 

BCS 

0.000 

2.20 

2.05 

2.50 

0.100 

2.19 

2.04 

2.49 

0.150 

2.17 

2.02 

2.48 

0.200 

2.14 

2.00 

2.46 

0.270 

2.09 

1.95 

2.41 

0.320 

2.04 

1.90 

2.37 

0.420 

1.91 

1.74 

2.25 

0.520 

1.72 

1.51 

2.10 

0.620 

1.40 

1.15 

1.86 

0.720 

0.86 

0.45 

1.49 

0.750 

0.62 

0.00 

1.34 

0.780 

0.30 

- 

1.15 

0.800 

0.00 

- 

1.00 

0.850 

- 

- 

0.55 

0.900 

- 

- 

0.00 


(iv) Effect of small external field (h = 0.4 T) is also studied on superconductivity. 
We observe from figs. 1 and 2 that small external field suppresses 
superconductivity. 

4.2 Magnetization 

Expression (29) is the expression from which we can study the variation of 
magnetization (m) with temperature for ferromagnetic superconductors. Using 
the values of various parameters given in tables 1 and 2 for UCoGe and Ulr 
and solving numerically, we obtain values of m at different temperature as 
given in tables 5 and 6 and the variation of m with T as shown in figs. 3 and 4 
respectively. From these graphs, we note that 

(i) Curie temperature for UCoGe and Ulr is 3K and 46 K respectively. 

In addition to the solutions described above, we have solved equations 
(26) and (34) for the coexistent state (A ^ 0,1 ^ 0). Values of A and I within 
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Table 4: Superconducting order parameter (A) for Ulr 


Temperature (T) 

K 

Superconducting order parameter 

A = x x 10- 21 J 

Theory 

h = 0.4 T 

BCS 

0.000 

1.400 

1.250 

1.700 

0.010 

1.400 

1.240 

1.700 

0.020 

1.400 

1.230 

1.690 

0.030 

1.380 

1.200 

1.676 

0.040 

1.350 

1.160 

1.651 

0.050 

1.301 

1.094 

1.614 

0.060 

1.247 

1.018 

1.553 

0.070 

1.160 

0.924 

1.486 

0.080 

1.072 

0.806 

1.393 

0.090 

0.960 

0.676 

1.300 

0.100 

0.819 

0.520 

1.180 

0.110 

0.640 

0.350 

1.050 

0.120 

0.452 

0.165 

0.880 

0.130 

0.235 

0.000 

0.680 

0.140 

0.000 

- 

0.470 

0.150 

- 

- 

0.230 

0.160 

- 

- 

0.000 
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Figure 1: Variation of superconducting order parameter (A) with temperature (T) for 
UCoGe. For 1 = 0, BCS curve is shown. 
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Figure 2: Variation of superconducting order parameter (A) with temperature (T) for 
Ulr. For 1 = 0, and h = 0, BCS curve is shown. 


Table 5: Magnetization (m) for UCoGe 


Temperature (T) 
K 


Magnetization m (T) 


0.000 

0.300 

0.500 

0.700 

1.000 

1.200 

1.500 

1.700 

2.000 

2.300 

2.600 

2.800 

3.000 


1.9900 

1.9740 

1.9550 

1.9400 

1.9000 

1.8600 

1.7700 

1.6920 

1.5580 

1.3530 

0.9740 

0.5960 

0.0000 


the coexistent state for UCoGe are given in table 7 and behaviour of A and 
I is shown in fig. 5. The behaviour of A and I in the coexistence region for 
UCoGe is quite unusual. 

Analysis of the expression (29) shows that the presence of magnetization 
reduces the volume of the phase space that is available for the Cooper 
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Table 6: Magnetization (m) for Ulr 


Temperature (T) 

K 

Magnetic order parameter 

m (T) 

0.0 

1.750 

5.0 

1.740 

10.5 

1.720 

16.8 

1.650 

24.1 

1.520 

29.5 

1.349 

33.2 

1.200 

37.0 

1.000 

40.5 

0.700 

45.0 

0.190 

46.0 

0.000 
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Figure 3: Variation of Magnetization with Temperature for UCoGe. 


pair. This leads to a decrease in the superconducting gap (A) and hence T. 
also reduces compared to the BCS value for given values of other relevant 
parameters. Obviously, increase in magnetization pushes the gap away from 
the Fermi surface and reduces its magnitude also. Our results agree with 
Karchev et al. [56] and Dahal et al. [57]. 
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Figure 4: Variation of Magnetization (m) with Temperature for Ulr. 

Table 7: Order parameters (A, I) in the ferromagnetic, superconducting and 
ferromagnetic - superconducting coexistent states for UCoGe. 


Temperature (T) 

K 

A 

coex 

1 

coex 

I 

A 

0.00 

1.90 

2.20 

0.40 

0.25 

0.20 

1.89 

2.15 

0.46 

0.16 

0.40 

1.87 

1.90 

0.50 

0.12 

0.60 

1.85 

1.13 

0.40 

0.16 

0.80 

1.80 

0.24 

0.00 

0.26 

0.85 

1.78 

0.00 

- 

- 

1.00 

1.72 

- 

- 

- 

1.20 

1.63 

- 

- 

- 

1.40 

1.52 

- 

- 

- 

1.60 

1.40 

- 

- 

- 

1.80 

1.24 

- 

- 

- 

2.00 

1.05 

- 

- 

- 

2.20 

0.85 

- 

- 

- 

2.40 

0.65 

- 

- 

- 

2.60 

0.45 

- 

- 

- 

3.00 

0.00 

- 

- 

- 


56 















Theoretical Study 
of Interplay 
Between 
S uperc onduc tivity 
and Itinerant 
Ferromagnetismy 


Figure 5: Order parameters of the ferromagnetic state (A = 0,1), the superconducting 
state (A , I = 0) and the coexistence state (A ^ 0,1 ^ 0) as a function of temperature 
for UCoGe. 



Figure 6: The energy spectrum of ferromagnetic superconductors. The upper curve 
corresponds to equation (30) and the lower one for equation (31). The upper horizontal 
line is at Um + A and the lower horizontal line is at Um - A. The gap is 2A. 


4.3 ENERGY SPECTRA AND DENSITY OF STATES 

If we compare the quasi particle energy spectra (Equations (30) and (31)) with 
the BCS energy spectrum, E = +A 2 , we sees that the modification is 

due to the presence of ferromagnetic energy. The plot of energy spectrum is 
shown in fig. 6. It is very clear that the gap is not at the Fermi level, rather it is 
pushed up. Now, the energy spectrum is symmetric around Um, contrary to the 
BCS case where it was symmetrical around E = 0. 
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to (j: x 10~ 31 J) 

Figure 7: Density of States for UCoGe for A = 1.15 x 10 2I J and A = 1.50 x 10 21 J. 



Figure 8: Variation of Density of states with u for Ulr for A = 1.10 x 10 21 J and A = 
1.30 x 10 2I J. 


Expressions (40) and (41) converge to the density of states of the BCS model 
for I = 0 and h = 0, i.e. 



(for |uj| > A) 
(for |uj| < 0) 


Solving numerically equation (40) for UCoGe and Ulr, we obtain the variation 
of N(to) / N(o) with to as shown in figures 7 and 8 respectively and corresponding 
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Table 8: Density of states for UCoGe 


u 



N (lo)/N(0) 


(x x 10- 21 J) 

BCS 

A = 1.15 x 10- 21 J A = 1.50 x 10 21 J 

h = 0.4 T 

0.10 

- 

2.75 

- 

- 

0.20 

1.400 

2.42 

2.85 

2.99 

0.40 

1.087 

1.90 

2.25 

2.55 

0.60 

1.000 

1.55 

1.85 

2.14 

0.80 

1.000 

1.34 

1.60 

1.45 

1.82 

1.62 

1.00 

1.000 

1.24 

1.30 

1.45 

1.50 

1.000 

1.17 




1.25 

1.37 

2.00 

1.000 

1.15 

1.25 

1.37 

2.50 

1.000 

1.15 

1.25 

1.38 

3.00 

1.000 

1.15 
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Table 9: Density of states for Ulr 


10 


N (u)/N(0) 


(x x 10- 21 J) 

BCS 

A = 1.10 x 10 21 J 

A = 1.30 x 10 21 J 

0.30 

1.40 

2.80 

3.50 

0.40 

1.20 

2.03 

2.60 

0.60 

1.05 

1.55 

1.85 

0.80 

1.00 

1.33 

1.60 

1.00 

1.00 

1.25 

1.45 

1.50 

1.00 

1.18 

1.35 

2.00 

1.00 

1.16 

1.30 

2.50 

1.00 

1.15 

1.28 

3.00 

1.00 

1.15 

1.27 


values are recorded in tables 8 and 9. We note that gap is pushed up (compared 
to BCS gap), giving rise to finite density at the Fermi level. Looking at the 
energy spectrum and density of states, one finds that the physical properties of 
the system will be dominated by the normal metal like behaviour. Our results 
agree with Dahal et al. [57], and Karchev et al. [56]. Small external field also 
pushes up A and in turn density of states rise at the Fermi level. 

4.4 SPECIFIC HEAT 

Solving numerically the expression (37) for specific heat divided by temperature 
(C/T) for UCoGe and Ulr, values obtained are given in tables 10 and 11 and 


59 













Kakani, S. 
Nuwal, A. 
Kakani, S.L. 


Table 10: Specific heat divided by temperature (C/T) for UCoGe 


Temperature (T) 

K 

C/T (mj / mol 

-K 2 ) 

BCS 

Theory Experimental 

0.10 

0.90 

10 

40.00 

0.20 

2.00 

20 

45.71 

0.30 

5.00 

30 

51.43 

0.40 

8.50 

40 

57.14 

0.50 

13.00 

50 

62.86 

0.60 

20.00 

60 

68.57 

0.70 

30.00 

70 

74.29 

0.80 

45.00 

80 

80.00 

1.00 

80.00 

- 

- 


Table 11: Specific heat divided by temperature (C/T) for Ulr 


Temperature (T) C/T (mj / mol - K 2 ) 

K 

Theory 

BCS 

0.02 

8.45 

0.90 

0.04 

16.90 

2.50 

0.06 

25.35 

5.20 

0.08 

33.80 

10.00 

0.10 

42.25 

16.50 

0.12 

50.70 

26.00 

0.14 

60.00 

40.20 

0.16 

- 

60.00 


variation of C/T versus T is obtained as shown in figures 9 and 10 respectively. 

From the graphs, we note: 

(i) A linear temperature dependence of C/T in coexistent state of 
ferromagnetism and superconductivity obtained from theory is in excellent 
agreement with experimental results for UCoGe [58] as opposed to the 
exponential decrease of specific heat in the BCS theory. Our results agree 
with Karchev et al. [56]. 

(ii) Specific heat studies also confirm that T c of ferromagnetic superconductors 
is affected by magnetic field, i.e. decreases slightly as compared to BCS 
value. 


60 














Theoretical Study 
of Interplay 
Between 
S uperc onduc tivity 
and Itinerant 
Ferromagnetismy 


Figure 9: Variation of Electronic Specific heat divided by temperature (C/T) for 
UCoGe. Experimental curve [4] . 



Figure 10: Variation of Electronic specific heat divided by temperature (C/T) for Ulr. 

4.5 FREE ENERGY AND CRITICAL FIELD 

From equation (55), we note that 

F fs “ F nf * 0 and A2 * 0 

This clearly reveals that the phase transition from normal ferromagnetic 
state to superconducting -ferromagnetic state is first order phase transition 
[58], 
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Table 12: Free energy difference (F sp - F nf ) in the coexistent state for UCoGe 


T/T c Free energy difference 
- (F - F ) 

_ sr _ 1 nt' _ 


0.00 

0.480 

0.10 

0.470 

0.20 

0.443 

0.30 

0.400 

0.40 

0.350 

0.50 

0.260 

0.60 

0.139 

0.70 

0.064 

0.80 

0.023 

0.90 

0.005 

1.00 

0.000 


Table 13: Free energy difference (F sf - F nf ) in the coexistent state for Ulr 

T/T. Free energy difference 

-< F s,--F nf > 


0.00 

0.45 

0.10 

0.44 

0.20 

0.41 

0.30 

0.37 

0.40 

0.32 

0.50 

0.23 

0.60 

0.10 

0.70 

0.04 

0.80 

0.02 

0.90 

0.01 

1.00 

0.00 


Expression (53) also shows that superconducting - ferromagnetic (SF) 
state has lower energy than the normal ferromagnetic (NF) state and therefore 
realized at low enough temperature. Our results agree with Karchev [56]. 
The free energy difference as function of temperature is shown in tables 
12 and 13 and variation shown in figures 11 and 12 for UCoGe and Ulr 
respectively. From these graphs, we note that the coexistent phase of 
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Figure 11: Free energy difference (F sp - F nf ) as function of T/T c for UCoGe, which 
displays coexistence of ferromagnetism and superconductivity below T = 0.8 K. 



T/T t 


Figure 12: Free energy difference (F sf - F nf ) as function of T/T c for Ulr, which displays 
coexistence of ferromagnetism and superconductivity below T c = 0.14 K. 

superconductivity and ferromagnetism is energetically favored compared to 
normal ferromagnetic case, which is consistent with the experimental fact 
that a transition to superconductivity - ferromagnetic state occurs well below 
the Curie temperature [4,59]. 
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Table 14: Critical Magnetic Field (H c ) for UCoGe 


Temperature (T) 

Critical Magnetic Field 

K 

(H (T)) 

0.00 

1.600 

0.10 

1.590 

0.20 

1.523 

0.30 

1.413 

0.40 

1.232 

0.50 

0.950 

0.60 

0.616 

0.70 

0.300 

0.80 

0.000 

Table 15: Critical Magnetic Field for Ulr 

Temperature (T) 

Critical Magnetic Field 

K 

H(T) 

0.00 

1.00 

0.02 

0.97 

0.04 

0.92 

0.06 

0.84 

0.08 

0.74 

0.10 

0.60 

0.12 

0.39 

0.14 

0.00 


Variation of critical field with temperature obtained numerically by solving 
equation (57) is shown in table 14 and 15 for the systems UCoGe and Ulr 
and variation is shown in figures 10 and 11 respectively. Critical field curves 
reveal that physical properties of the system are dominated by normal metal 
like behaviour. 


5. CONCLUSIONS 

In the BCS theory of superconductivity, the conduction electrons in a 
metal cannot be both ferromagnetic ally ordered and superconducting. 
Superconductors expel magnetic fields passing through them but strong 
magnetic fields kill the superconductivity. Even small amounts of magnetic 
impurities are usually enough to eliminate superconductivity. Much work has 
been done both theoretically and experimentally to understand the interplay 
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between ferromagnetism and superconductivity and to search for the possibility 
of true coexistence of these two ordered states. Prior to 2001, in all known 
ferromagnetic superconductors, the superconductivity phase is observed in a 
small part of the phase diagram. This situation has changed with the emergence 
of superconductivity under the background of ferromagnetic state in ZrZn 2 , 
UGe 2 , URhGe, UCoGe and Ulr. In these ferromagnetic superconductors, 
a superconducting transition takes place at a temperature T c deep in the 
ferromagnetic state, i.e. well below the Curie temperature T pM , without 
expelling magnetic order. In uranium intermetallic compounds magnetism 
has a strong itinerant character and both ordering phenomena are carried by 
5f electrons. Thus, it is proper to study a model where the coexistence of 
both ferromagnetism and superconductivity can be described by only one 
kind of electrons. Such a model study has recently been initiated by Karchev 
et al. [56]. 

We have studied the inteiplay between itinerant ferromagnetism and 
superconductivity in a model single - band homogeneous system by using 
a mean - field approximation. The superconducting part is treated by BCS 
theory and the itinerant ferromagnetic part is studied with the use of Hubbard 
Hamiltonian. Following Green’s functions technique and equation of motion 
method, we have obtained coupled equations of superconducting gap (A) 
and magnetization (m). The two order parameters, m and A are coupled with 
each other through the distribution functions, which are also functions of m 
and A. We found that there are three nontrivial sets of solutions for A and I 
or m : the ferromagnetic (A = 0, m ^ 0), the superconducting (m = 0, A ^ 
0), and the coexistent state (A ^ and m = 0). 

We have closely studied the specific heat, energy spectra and density of 
states and free energy for a ferromagnetic superconductor considering that 
same electrons are responsible for both ordered states and superconductivity 
is due to s - wave pairing. 

The specific heat has a linear temperature dependence (C oc T) in the 
coexistence region as opposed to the exponential decrease of the specific 
heat in the BCS theory. Our result fits well the data observed experimentally 
for UCoGe. The specific heat has a linear temperature dependence as in 
the normal ferromagnetic metals, but increases anomalously at small 
magnetizations. 

Free energy studies shows that the superconducting-ferromagnetic 
coexistence state has lower energy than normal ferromagnetic state and 
therefore coexistence state is realized at low enough temperature. Moreover, 
the ferromagnetic state never appears when the superconducting state is 
already in existence. This behaviour seems to be a particular feature of 
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Figure 13: Variation of Critical Magnetic Field (H c ) with Temperature for UCoGe. 


the itinerant -electron model. For T c < T FM , there is a ferromagnetic to 
superconducting-ferromagnetic transition at T = T < T , where T is the 
temperature below which superconductivity and ferromagnetism coexist. 
For T pM < T. and T < T , the superconducting state is more stable than the 
ferromagnetic state. This means that the ferromagnetic state can never be 
reentrant when the superconducting state is already present [10,60-63]. This 
is contrary to the situation in rare earth ternary compounds ErRh 4 B 4 and 
HoMo 6 S g [60,64], where superconductivity occurs at higher temperatures, 
and at lower temperatures the ferromagnetic state may appeal'. 

This may be physically understood as follows: for localized ferromagnetism, 
the saturation magnetization, which depends on the magnitude of an 
individual local moment, and the Curie temperature, which depends on the 
exchange interaction between neighboring moments are not related to each 
other. The energy reduction due to the saturation magnetization may be quite 
large while the value of T FM can be very small. For an itinerant system, the 
situation is quite different. Both the saturation magnetization and T pM are 
decided dominantly by the Hubbard interaction U and the band structure. 
For a system with low T pM , both the saturation magnetization and the energy 
reduction due to the magnetization will certainly be small. Therefore, when, 
T fm < T , the itinerant ferromagnetic state is always at a disadvantage in 
energy as compared with the superconducting state. 
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Figure 14: Variation of critical magnetic field (H) with Temperature (T) for Ulr. 



Figure 15: Phase diagram of itinerant ferromagnetic superconductors. 

Blount and Varma [64] pointed that the electromagnetic coupling which 
results from the interaction between the superconducting d electrons and 
local magnetic moments of rare - earth ions plays a crucial role in ternary rare 
earth superconductors, ErRh 4 B 4 and HoMo 6 S g . For itinerant model, however, 
both the magnetism and superconductivity come from the electrons in the 
same f band. Here the main coupling is the Coulomb interaction between 
electrons, and the Hubbard model is an approximate description of this 
interaction. In our work, we have treated the magnetization as uniform and 
thus the electromagnetic coupling is absent. 
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The study of the energy spectrum and the density of states reveals that the 
physical properties of the system are dominated by normal metal l ik e behaviour. 

The itinerant electron model developed by us provides a microscopic 
relationship between ferromagnetism and superconductivity. This model 
explains satisfactorily the behaviour of specific heat, energy spectra and 
density of states, free energy, etc. Our study unambiguously shows that 
superconductivity and itinerant ferromagnetism truly coexist in these systems 
and both ordering phenomena are carried by the same 5f electrons. Free 
energy study clearly reveals that it is possible to become superconducting 
via a first - order phase transition if the system on cooling first shows 
ferromagnetism. Our results are in agreement with Karchev [55], Linder and 
Sudbo [54] and Dahal et al. [57]. 

A possible phase diagram based on our model is as shown in fig. 15. 

It is well known that mean - field approximation has its limitation in 
treating ferromagnetism. We are not sure, whether and how an electromagnetic 
coupling or the coupling between superconducting electrons and magnetic 
fluctuations will be generated if the relevant part of the Coulomb interaction 
between f electrons is correctly taken into account. This question is certainly 
important and needs further investigation. For a more rigorous approach, the 
effect of magnetic fluctuations should also be considered. Since the mean 

- field theory is qualitatively correct away from T , we expect that the 
essential features of our results should remain valid. 

We hope in the near future measurements of the electronic and magnetic 
excitation spectra in the superconducting and magnetic phases of these 
systems will reveal crucial information on the superconducting gap structure 
and pairing mechanism. UCoGe may be the first material to reveal proof for 
the existence of the long - searched spontaneous vortex phase. 

Progress in this field requires a new generation of pure crystals. Up to now, 
the attempts to discover Ce ferromagnetic superconductivity have failed. 

With the discovery of weak itinerant ferromagnetic superconductors a 
new research theme in the field of magnetism and superconductivity has 
been disclosed. Research into ferromagnetic superconductors will help to 
unravel how magnetic fluctuations can stimulate superconductivity, which is 
a central theme running through materials families as diverse as the heavy 
fermion superconductors, high - T c cuprates and recently discovered FeAs 

- based superconductors [65,66]. This novel insight might turn out to be 
crucial in the design of new superconducting materials. 

An interesting aspect of ferromagnetic superconductivity concerns the 
influence of ferromagnetism on macroscopic phenomena such as Meissner 
effect [38]. Other related topics for further study are superconductivity in 
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ferromagnetic domain walls and the phenomena associated with the relative 

orientation of the superconducting order parameter to the magnetization, the 

effect of superconductivity on ferromagnetic domain structure [67], etc. 
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